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Abstract 

The  Justification  of  continuing  the  Rayleigh  expansion  to  the  grating's 
surface  (the  Rayleigh  hypothesis)  and  its  convergence  properties  are  considered. 
A  class  of  gratings  for  which  the  Rayleigh  hypothesis  is  exact  is  identified,  a 
prime  example  'of  which  is  the  sinusoidal  grating  (SG) .  Based  on  an 
identification  of  the  origin  of  the  limited  stability  of  the  Rayleigh  expansion, 
a  modified  expansion  is  introduced,  dubbed  as  the  dressed  Rayleigh  expansion. 
This  new  expansion  presumably  has  excellent  convergence  properties  as  explicitly 
demonstrated  for  the  sinusoidal  grating.  The  dimensionality  N  of  the  matrix 
which  must  be  inverted  for  a  SG  of  arbitrary  depth  g  and  periodicity  d  is  found 
to  be  N  -  8vg/d. 
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1 .  Introduction 

Light  scattering  from  a  grating  is  hardly  a  new  phenomenon,  and  consequently 

1  2 

the  literature  on  the  subject  is  enormous.  ’  It  is  important  however  to 
realize  that,  barring  few  recent  exceptions,  the  majority  of  the  literature 
deals  with  very  shallow  gratings,  i.e.,  where  the  ratio  between  the  height  g  and 
periodicity  d  is  g/d  <<  1.  Shallow  gratings  are  realized  in  many  important 

physical  situations,  such  as  optical  refraction  gratings,  holograms,  gratings 

3  4 

induced  py  surface  phonons,  and  a  large  class  of  volume-gratings. 

Studies  of  light  scattering  from  deep  gratings,  i.e.,  when  g/d  «  1  or 

5-1 1 

larger,  are  relatively  recent.  The  motivation  to  consider  this  regime  stems 

in  part  from  the  advent  of  new  fabrication  capabilities.  By  using  holographic 

12  13 

techniques,  or  laser-induced  deposition  from  a  volatile  organo-metallic  gas, 

it  is  possible  to  form  gratings  with  g/d  £  1 .  Another  reason  for  the  interest 

in  deep  gratings  is  the  expected  new  qualitative  features.  When  B  «  1,  where 

hereafter  we  use  0  »  2vg/d  as  the  natural  parameter,  the  surface  can  only  to  a 

limited  degree  exchange  surface-parallel  momentum  quanta  kG  -  2ir/d  with  the 

incident  light.  Consequently,  for  shallow  gratings  it  is  sufficient  to 

consider  one  or  two  Bragg  reflections  in  addition  to  the  dominating  specularly 

reflected  and  transmitted  waves.  For  deep  gratings  (B  >>  1),  the  situation  is 

quite  different.  Here  the  grating  can  efficiently  exchange  with  the  incident 

light  a  large  number  of  quanta  k^.  Therefore,  the  reflected  wave,  for  example, 

is  made  up  of  many  Bragg  reflections,  strongly  interfering  with  each  other  and 

with  the  specularly  reflected  wave.  Thus,  unlike  the  0  <<  1  regime  which  can  be 

described  in  terms  of  few  weak  Bragg  reflections  superimposed  on  a  strong 

specularly  reflected  wave,  when  0  >>  1  the  identity  of  these  components  is 

blurred  and  the  reflected  light  acquires  a  new  character. 
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The  qualitative  distinction  between  the  g  <<  1  and  g  >>  1  regimes  is 

reminiscent  of  the  situation  with  regard  to  band  structure  in  Condensed  Matter 
1 4 

Physics.  The  umklapp  scatterings  play  the  role  of  Bragg  reflections,  and  the 
regimes  of  weak  and  strong  periodic  ionic  potentials  are  the  counterparts  of  the 
shallow  and  deep  grating  regimes,  respectively.  In  keeping  with  this  analogy, 
the  grating's  surface  modes  (surface  plasmons)  correspond  to  the  band  (bound) 
states.  By  the  same  token,  we  can  draw  an  analogy  between  light  scattering  from 
a  deep  grating  vs.  from  a  rough  surface  to  the  electronic  states  in  ordered  and 
disordered  systems,  respectively.  In  both  cases,  the  difference  is  the 


underlying  symmetry  of  the  system:  a  discrete  translational  symmetry  in  gratings 
.(crystals)  and  none  at  all  in  rough  surfaces  (disordered  systems).  As  a  result, 
the  character  of  the  ensuing  interference  among  the  many  reflected  components  is 


sri 


fundamentally  different  in  both  cases. 

Attempts  to  study  theoretically  the  deep  grating  regime  have  encountered 

1  15 

numerical  convergence  problems  V  when  g/d  >  0.1.  The  associated 
1  2  1 6  IT 

literature  ’  ’  V  often  makes  reference  to  two  related  issues  which  are  at  the 

focus  of  this  work:  the  validity  of  the  Rayleigh  hypothesis  and  the  convergence 

1 8 

of  the  Rayleigh  expansion.  The  Rayleigh  expansion  of  the  electromagnetic 
fields  above  the  selvedge  domain  (region  (a)  in  Fig.  I.a)  embodies  the 

1  M 

underlying  symmetry  of  the  grating  as  expressed  by  the  Floquet-Bloch  theorem, 
and  the  boundary  conditions  for  z  •*  •  (in  the  notation  of  Fig.  I.a).  Thus  the 
Rayleigh  expansion  results  from  a  general  principle.  The  first  issue  is 
whether,  or  when,  it  can  be  continued  into  the  selvedge  domain  all  the  way  to 
the  grating's  surface  (domain  (b)  in  Fig.  I.a).  The  assumption  that  this  is 
possible  is  referred  to  as  the  Rayleigh  hypothesis.  The  second  issue  is  what 
are  the  convergence  properties  of  the  Rayleigh  expansion,  in  light  of  the 
numerical  difficulties  encountered  in  applying  it. 


Notwithstanding  the  many  attempts  to  settle  these  issues  for  a  general 
grating,^ »2»l6»17t19  there  n0  ciear  answer  yet.  This  state  of  affairs  is 
particularly  puzzling  in  view  of  the  generality  of  the  Rayleigh  expansion  (which 
is  actually  a  Fourier  series  over  a  unit  cell  of  length  d) .  In  an  effort  to 
understand  these  issues,  we  focus  on  the  analysis  of  light  scattering  from  a 
class  of  simple  gratings,  with  the  sinusoidal  grating  (SG)  serving  as  a 
prototype. 

The  two  results  of  this  work  are  the  following: 

(i)  We  identify  a  class  of  gratings  for  which  the  Rayleigh  hypothesis 
is  exact.  Gratings  in  this  class,  such  as  the  SG,  have  the  property  that  the 
domain  above  the  selvedge  and  the  adjacent  domain  in  the  selvedge  (e.g. ,  domains 
"a"  and  "b"  in  Fig.  I.a),  share  the  whole  x-axis  except  for  a  set  of  isolated 
nonsingular  points. 

(ii)  With  regard  to  the  convergence  properties  of  the  Rayleigh 
expansion  convergence,  we  identify  the  cause  for  divergence  when  3  >>  1.  This 
analysis,  in  turn,  leads  naturally  to  a  simple  alternative  expansion  —  the 
dressed  Rayleigh  expansion  of  Eqs.  (3-5)— (3-6)  —  which  has  presumably  excellent 
convergence  properties.  This  proposition  is  explicitly  checked  for  the  SG.  We 
find  in  this  case  that  the  dressed  expansion  converges  for  an  arbitrary  3,  and 
the  estimated  ,order  of  the  matrix  needed  to  be  inverted  Is  M3.  Hence,  for 
example,  for  a  SG  with  g/d  -  1,  inverting  a  matrix  of  the  order  £  50  should  be 
adequate. 

The  paper  is  organized  as  follows:  In  Section  2  we  introduce  the  Rayleigh 
expansion  and  discuss  the  class  of  gratings  for  which  the  Rayleigh  hypothesis  is 
exact.  Section  3  is  devoted  to  examining  the  convergence  properties  of  the 
Rayleigh  expansion  and  the  introduction  of  the  dressed  expansion.  In  Section  M 


y.  v 
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the  latter  is  analyzed  for  the  SG.  Discussion  and  concluding  remarks  are  given 
in  Section  5. 

2.  Rayleigh  Expansion  and  Hypothesis 

To  introduce  the  Rayleigh  expansion,  we  start  by  considering  the  symmetry  of 

the  system:  Since  the  grating  is  invariant  under  translations  of  the  type 

x  ■*  x  +  Ld,  where  L  is  any  integer  and  d  is  t  •  periodicity  (all  notations  are 

defined  in  Fig.  I.a),  so  are  the  solution  of  Maxwell's  equations.  Consequently* 

6  B  i  ]  y 

according  to  the  Floquet-Bloch  theorem,  *  ’  *  the  field  eigenmodes 

( electric/raagnetic  fields)  1*(x,z)  satisfy 

ik„d 

’F(x-*-d,z)  -  e  ¥(x,z)  ,  (2.1) 

where  k  is  a  surface  parallel  (in  the  x-direction)  momentum  label.  For  a 

scattering  wave  mode,  k^  is  the  x-eomponent  of  the  incident  wave.  For  a  surface 

wave  mode  (surface  plasraons  in  the  case  of  metallic  gratings),  k|(  is  a 

continuous  label  restricted  to  the  first  Brillouin  zone  -kG/2  S  k(J  S  kQ/2*  To 

simplify  the  subsequent  analysis,  we  consider  only  the  configuration  of  a  p- 

polarized  plane  wave  incident  perpendicular  to  the  grating's  grooves  direction 

(the  y-axis) .  In  this  case,  a  diagram  in  the  xz-plane  is  sufficient. 

The  general  symmetry  property  (2.1)  determines  the  functional  form  of  the 

fields  up  to  constants,  to  be  determined  by  matching  the  proper  boundary 

"1“ot 

conditions.  Assuming  an  e  time  dependence  of  the  electromagnetic  fields, 

where  Wq  is  the  field  frequency,  we  can  always  expand  the  solutions  of  Maxwell's 

equations  (at  least  when  they  are  not  singular)  in  terms  of  the  following 

8  1 6  1 8 

complete  set  of  functions  in  the  0  <  x  <  d  interval:  *  * 

•  ^  i[k.x  -  W  U)z]  .  i[k.x  +  W  (i)z] 

2  (X.z)  -  l  lc  (l)p  _(l)e  *■  a  +  A  (l)p  (l)e  a  } 

U  £a.(e  U  U 


'  * 


(2. 2. a) 


§  (x,z)  «  s  l  {c  (£)< 


iLk  x  -  W  U)z] 

X  Ct 


+  A  (i)e 
a 


i[V  «■  WoU)z]|  f 


where 


( 2.2. b) 


k(o)  *  /I~  k  , 
a 


k  *  Wq/c  , 


k.  -  2tt/ d 
u 


wa(i)  -  [k2(a)-k2]1/2  .  k£  -  k  +  lkQ 


A  A  A  A 


>.,±  ‘  -kuT  Ckiz  ♦  V11*1  •  s '  *  "  2 


(2.2.C) 


In  (2.2)  we  choose  always  IraCfcMfc)]  or  ReCW^d)]  as  positive,  x,  z  and  s  are 

unit  vectors  in  the  x-,  z-  and  y-directions  (Fig.  I.a),  and  a  denotes  a  domain 

in  the  xz-plane  with  a  constant  dielectric  function  e  .  (e  can  depend  on 

a  a  0 

in  the  following  analysis  since  is  kept  fixed  throughout.)  Once  a  convenient 
partition  of  the  xz-plane  into  domains  has  been  chosen,  the  coefficients  C^.A^ 
are  determined  by  matching  the  boundary  conditions  across  the  boundaries  of  the 
domains. 

To  demonstrate  the  choice  of  domains  a,  we  consider  in  particular  the  SG- 
type  example  (Fig.  I.a).  It  is  obvious  that  in  this  case  there  are  four 
relevant  domains  and  hence  four  expansions  of  the  type  (2.2)  to  determine.  In 
these  terms,  the  issue  of  whether  the  Rayleigh  expansion  can,  or  cannot,  be 
continued  to  the  grating’s  surface  is  tantamount  to  whether  the  expansions  in  domains 


’a”  and  "b' 


(and  similarly  "c"  and  "d")  are  identical.  This  issue  must 


be  decided  by  matching  the  boundary  conditions  between  domains  "a"  and  "b"  along 
the  separating  plane  z-g. 

Consider  for  example  the  tangential  components  of  £  and  B  along 

A  A 

the  x  and  s  directions  [Eq. (2.2)].  By  comparing  coefficients  of  the 


ikix, 

complete  set  (e  }  we  obtain 

C  U)e"R  -  A  U)e+R  -  C.  U)e~R  -  A.  U)e+R 
a  a  b  b 

C  U)e"R  +  A  U)e+R  =  C.  U)e“R+  A.  U)e+R  (2.3) 

a  a  b  b 

where 

R  =  iW  (£)g  -  iW.  U)g  .  (2.4) 

3  D 

Equations  (2.3)  yield  C  (li.)  -  C^d)  and  Aa(i,)  -  A^ ( il ) ,  and  hence  the  Rayleigh 
hypothesis  is  exact.  In  establishing  (2.3)  we  used  the  fact  that  domains  "a" 
and  "b"  share  the  entire  x-axis  (or  unit  cell).  Thus  it  is  possible  to  equate 

ikgX 

coefficients  of  each  member  of  the  complete  set  {e  }.  It  has  been  tacitly 

assumed  at  this  junction  that  the  set  of  isolated  nonsingular  points  (x  «  nd, 

z  -  g),  jointly  shared  by  domains  "a",  "b"  and  "c”,  do  not  break  the  orthogonality  of 
ik  x 

the  set  {e  *  }  over  a  unit  cell.  The  above  reasoning  remains  valid  for  any 
grating  such  that  the  selvedge  domain  and  the  domain  above  it  (e.g. ,  domains  "a" 
and  "b"  in  Fig.  I.a)  share  the  entire  x-axis  except  for  a  set  of  isolated 
nonsingular  points.  Figure  l.b  depicts  a  slightly  more  complicated  example 
where  the  Rayleigh  expansions  in  domains  "a"  and  "b"  are  identical. 

An  example  for  which  the  above  argument  fails  is  the  square-well  grating, 

Fig.  I.c.  Here  again  we  start  by  dividing  the  xz-plane  into  four 
domains  and  try  to  match  domains  "a"  and  "b".  It  is  obvious  that  since  the 
boundary  between  domains  "a"  and  "b"  is  not  the  entire  x-axis,  Eq.  (2.3)  is  not 
valid;  hence  the  Rayleigh  hypothesis  is  not  exact.  For  a  shallow  square-well 
•grating,  however,  the  Rayleigh  hypothesis  may  provide  a  good  approximation. 

The  foregoing  discussion  is  in  keeping  with  the  accepted  point  of  view  that 
in  general  the  Rayleigh  expansion  above  the  grating  cannot  be  continued  to  the 
grating’s  surface.  We  have  demonstrated,  nevertheless,  that  for  a  certain 
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class  of  gratings,  and  the  SG  in  particular,  the  Rayleigh  hypothesis  is  exact. 
Note  also  that  we  restricted  ourselves  to  nonsingular  fields  (and  first 
derivatives)  only.  When  singularities  are  present,  the  implication 
is  of  surface  charges,  an  infinite  amount  of  charge  at  particular  points, 

surface  currents,  etc.  Consequently,  the  above  conclusions,  which  are  based  on 
the  "no-surface-charges/currents"  boundary  conditions,  are  expected  to  be 
modified. 


3.  Dressed  Rayleigh  Expansion 

Even  when  the  Rayleigh  hypothesis  is  exact,  the  usefulness  of  the  Rayleigh 

expansion  must  be  substantiated  by  verifying  good  convergence  properties.  As 

1 5 

difficulties  in  applying  the  expansion  to  deep  gratings  indicate,  this  is 

probably  not  the  case.  It  is  for  this  reason  that  alternative  schemes  have  been 
1  2  6-8  1 6 

suggested  ’  ’  ’  with  a  presumably  wider  range  of  convergence.  Our  strategy, 

on  the  other  hand,  is  first  to  expose  the  deficiency  of  the  Rayleigh  expansior 
rather  than  to  abandon  it.  Once  this  is  achieved,  we  are  naturally  led  to  a 
modified  expansion  which  we  believe  (and  demonstrate  for  the  SG  in  the  next 
section)  to  have  excellent  convergence  properties  for  deep  gratings. 

When  the  Rayleigh  hypothesis  is  exact,  it  is  sufficient  to  consider  only  two 
domains  in  the  xz-plane  (see  Fig.  2).  Consequently  there  are  only  two 
expansions  to  deal  with  [Eq.  (2.2)],  indexed  as  "0"  and  "1".  The  electric  field 
expansions  are  therefore 


£Q(x,z) 


“  iCk.x  -  Wrt(Jl)z] 

I  {C  ( Jt)p  _U)e 
A—®  u 


~  i[  kpx  +  WnU)z] 

+  y*>Po,+U)e  \ 

.  “  i[k.x  -  W.U)z] 

Ej  (x,z)  -  l  C1  U)p1  _U)e 

Irn-m  ' 


(3.1. a) 


and  the  magnetic  field  expansions  are 


Bq(x,z) 


(x,z) 


-  ”  iCk5x  -  W_(£)z] 

s  I  |cQU)e  11  0 

£=-» 


i, / ,  \  -  “  i[k„x  -  W.  ( £)z] 

s  l  C  (£)e  *  1 

£=-« 


i[k  x  +  W  (£)z] 

+  AQ(£)e  } 

(3*1 .b) 


In  Eq.  (3.1) 


C0U) 


6£,0  E0p 


(3-2) 


where  Eg  is  the  arbitrary  amplitude  of  the  incident  p-polarized  light. 

Equation  (3.2)  embodies  the  boundary  condition  of  the  downward-propagating 
incident  wave  in  region  "0"  (Fig.  2),  and  (3.1)  incorporates  the  outgoing-wave 
boundary  conditions  at  |z|  •*  ®. 

The  deficiency  of  (3.1)  becomes  apparent  by  considering,  for  example,  the 
Ag(£)-term  in  (3 . 1 . t>) .  For  |£|  -*•  ®,  the  x-independent  factor  tend3  to 

iW  (£)z  _k  | £ jz 

A  (£)e  u  - ►  A  (£)e  u  (3.3) 

0  |£|+»  0 

where  (2.2. c)  was  used.  For  the  z  >  0  portion  of  domain  "0",  the  z-dependent 
factor  in  (3.3)  converges  exponentially  with  |£|.  However,  (3.1)  is  valid 
throughout  region  "0".  In  particular,  at  the  bottom  of  the  troughs  where 
-g  £  z  S  0,  the  z-dependent  factor  in  (3.3)  diverges  exponentially  with  |£|. 
Consequently,  in  order  to  keep  the  total  field  BQ(x,z)  finite  for  z  £  0,  it 
follows  that  the  exact  Aq(£)  must  converge  for  |£|  -*■  »  at  least  exponentially. 
[This  argument  fails  for  gratings  which  give  rise  to  singularities  in  the  fields 
or  in  their  derivatives.]  Therefore,  for  z  <  0  the  Rayleigh  expanion  of  EQ  and 
is  a  sum  of  many  terms,  most  of  which  (large  |£|)  are  products  of 
exponentially  large  numbers  (the  z-dependent  factors)  times  exponentially  small 
numbers  (the  Ag(£)'s).  However,  in  an  actual  calculation  there  are  always 


errors  in  the  calculated  AQU)'s.  These,  in  turn,  will  lead  to  large  variations 
in  the  total  fields,  since  the  AqOO's  are  multiplied  in  (3.1)  by  very  large 
exponents.  This  type  of  numerical  instability,  demonstrated  here  for  the  z  £  0 
portion  of  the  "0"  domain,  has  an  adverse  effect  on  other  portions  of  the  xz- 

1  9  i  f.  oi 

plane  by  virtue  of  the  extinction  theorem.  ’  ’  *  Therefore,  for  cases  where 

many  Bragg  reflections  (  | Jl |  * s )  contribute  to  (3-1),  i.e.,  for  deep  gratings,  the 
Rayleigh  expansion  is  intrinsically  unsuitable.  For  shallow  gratings,  where 
only  a  few  Bragg  reflections  contribute,  the  instability  just  described  does  not 
arise. 

The  above  deficiency  of  (3.1)  can  be  easily  remedied  in  the  following 
manner.  Consider  again,  for  example,  the  AQ(£)-term  in  (3 . 1 . b) .  It  can  be 
rewritten  as 

iW„U)z  iW_(  i.)  ( z+g) 

A0U)e  =  a0U)e  u  ,  (3.H) 

where  g  is  the  (positive)  minimum  of  the  grating  in  region  "0"  (see  Fig.  2)  and 
UqU)  =  A0U)exp[-iW0U)g].  Since  by  construction  z  +  g  £  0  throughout  region 
"0",  the  exponential  factor  on  the  RHS  of  (3.*0  never  diverges;  in  fact,  it 
always  converges  exponentially  (and  is  unity  at  the  grating’s  profile). 

Furthermore,  since  (3.1)  is  valid  throughout  region  "0",  including  at  the 
grating  surface,  the  exact  oiq ( & ) ' s  must  diminish  with  very  high  l  to  render  a  finite 
total  field  at  all  points.  Consequently,  the  relevant  exact  (^(JO's  are  not 
necessarily  exponentially  small,  and  hence  small  errors  in  the  calculated 
OgOD's  will  not  lead  to  instabilities  in  the  total  field.  These  considerations 
suggest  that  we  transcribe  the  Rayleigh  expansion  (3.1)  into  a  "dressed  Rayleigh 
expansion"  by  writing 
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iWnU)z  iW  (!)(z+g) 

AQ(4)g  -  aQ(£)e 


-iW.OOz  -1W  (!)(z-g) 

C  (l)e  1  **  Y,U)e  '  ,  (3.5) 


where  the  dressed  expansion  coefficients  are  given  by 

-iWQU)g 

a0U)  *  AQU)e 

-iv  (  9.W 

7,(11)  -  C,  (!)e  1  .  (3.6) 

The  central  assertion  of  this  work  is  that  the  dressed  Rayleigh  expansion, 
defined  in  (3.5)  and  (3.6),  has  excellent  convergence  properties  for  a  wide 
range  of  values  of  B  and  is  a  suitable  framework  for  deep  grating  calculations. 
This  proposition  is  explicitly  demonstrated  for  the  SG  in  the  next  section. 
Furthermore,  the  procedure  used  in  Section  4  for  the  SG  can  be  generalized  to 
other  gratings  (see  Section  5),  indicating  the  wide  applicability  of  the 

dressed  Rayleigh  expansion. 

4.  Sinusoidal  Grating  (SG) 

The  proposition  stated  as  a  conjecture  in  Section  3  with  regard  to  the 
dressed  Rayleigh  expansion,  (3.5)  and  (3.6),  is  now  explicitly  checked  for  the 
SG.  Our  starting  point  is  the  exact  infinite  set  of  coupled  linear  equations1^ 
satisfied  by  the  expansion  coefficients  of  (3-1)  (see  Appendix  A  for  the 
derivation) , 

]  "m.tV11  ‘  “B("> 

Jt*-" 


l  N l  .C. (!)  -  vB(m)  ,  (4.1) 

o— -  m* 1  1 


where  the  superscript  B  is  a  reminder  that  (4.1)  pertains  to  the  expansion 
coefficients  in  (3.1 )•  These  coefficients  are  hereafter  referred  to  as  "bare", 
in  distinction  from  the  "dressed"  expansion  coefficients  defined  in  (3.6). 
The  entries  to  (4.1)  are  explicitly  given  by  (see  Appendix  A): 


V-S<V0)  *  V">):lEo'p 


(1.2) 


where  the  symbol  J  denotes  the  Bessel  function  of  order  m,  and  all  other 
m 

R  R 

symbols  have  been  defined  in  (2.2).  The  Bessel  function  factor  in  M  and  N 

determines  the  mixing  between  the  Bragg  reflections  of  order  m  and  £,  and  hence 

involves  the  g-parameter.  Thus,  for  instance,  when  g  -  0  (flat  surface),  it 

follows  that  m  -  £,  i.e.,  no  mixing.  The  role  of  the  other  factors  in  (4.2)  is 
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elaborated  elsewhere. 

The  exact  equations  (4.1)  and  (4.2)  expose  from  yet  another  vantage  point 
the  deficiencies  of  the  Rayleigh  expansion  (3.1 ).  Consider  for  example  M  and 
the  corresponding  equations  in  (4.1).  When  m  is  fixed  and  |fc|  -*•  •,  the  fastest 
changing  factor  in  (4.2)  is  the  Bessel  function.  For  |£|  +  •  it  follows  that 
W^m)  <<  Wq(£)  -  ikQ|£|,  and  therefore,  up  to  uninteresting  factors  (phases  and 
powers  of  |£|),  the  matrix  elements  of  M  behave  as  J 

lim  M®  -  lira  J  U6|£j)  -  el*ln(6)  ,  (4.3) 

|£|-»  “•*  |  £  j  +•  1 

m  fixed 

where  again  6  -  k„g  and 

U 


n(B)  -  A  ♦  ln[ 8/(  1  +  /I'  *~P)]  . 


(4.4) 


Hence,  for  8  >>  1  the  matrix  elements  of  M  diverge  exponentially  with  |fc|.  On 


the  other  hand,  for  fixed  m  the  RHS  of  (4.1)  is  a  constant.  Consequently,  the 


exact  AqU)  must  converge  at  least  as  e  which  is  in  keeping  with  the 


analysis  of  Section  3-  Furthermore,  since  M  diverges  with  |S,|  as  indicated  in 
(4.3),  there  is  no  natural  point  at  which  to  truncate  (in  |H|)  the  matrix  and 


control  the  corrections  in  an  actual  calculation. 


The  situation  is  dramatically  changed  once  (4.1)  are  transcribed  to 


equations  for  the  dressed  coefficients,  Eq.  (3.6).  These  are 


“o(l)  ■ 


l  N°it  yt>  -  »%)  , 


(4.5) 


where  the  l-dressed  matrices  M°  and  ND  are 


D  B  iW0(1L)8 

m,  l  Mm,8.  6 


md  =,  mB  iW1(£)g 
m,  i.  m,  l  e 


(4.6) 


By  the  same  procedure  which  led  to  (4.3)  and  (4.4),  we  obtain 


lim  -  lim  J  [ig|i|]  e"8^!  -  e^^(e)  , 


m  fixed 


(4.7) 


where 


4»(8)  -  n(B)  -  B 


(4.8) 


and  we  have  ignored  in  (4.7)  uninteresting  factors  (phases  and  powers  of  (|l|). 
The  central  point  of  this  section  is  to  recognize  that  $(B)  is  negative  for  all 


values  of  8  (see  Fig.  3).  Moreover,  since 


me 


for  g  «  1 


♦  (B)  -  1  +  lng 
1_ 

~  "  28 


for  8  >>  1 


(4.9) 


we  have  an  estimate  for  the  dimension  N  of  the  effective  MD  matrix  (comprised  of 
the  significant  matrix  elements  in  MD).  For  8  >>  1  Eq.  (4.9)  yields 

N  *  2  *  28  =  8ng/d  ,  (4.10) 

where  the  extra  factor  of  2  in  (4.10)  is  added  to  account  for  positive  as  well 
as  negative  i-values.  The  estimated  N  in  (4.10)  is  roughly  the  number  of 
significant  terms  in  the  dressed  Rayleigh  expansion.  Phrased 

differently,  N  provides  the  natural  truncation  point  of  in  an  actual 
calculation.  Thus,  for  instance,  when  g/d  =  1  the  inversion  of  a  50  x  50  matrix 
should  be  quite  adequate. 

To  complete  the  analysis,  we  now  consider  how  the  convergence  of  (4.5)  for 
the  case  j  m J  ■*  <*>  and  51  is  determined.  Obviously,  the  m-convergence  factor  to  be  used 
is  arbitrary.  A  convenient  choice,  symmetric  to  (4.6),  is 


M  p  -  e 

fll  |  Kr 


i[W  (m)  +  W  (i)]g  iW  (m)g 

"5,4  ’  *(m>  '  e  U  (m) 


i[W  (m)  +  VLU)]g  R  iW  (m)g 

V*  “  e  \i  •  *■>  -  e  vB(m)  , 


(4.12) 


and  the  fully  dressed  equations  are 

CO 

J  V‘>  •  “(m) 

co 

l  N  Y^t)  -  v(n)  .  (4.13) 

The  set  (4.13)  is  exponentially  convergent  in  both  the  m-  and  ^-''directions" 
for  an  arbitrary  8.  It  is  therefore  a  perfectly  stable  framework  for 
calculations  of  deep  SO.  The  size  of  the  matrix  needed  to  be  inverted  is  of  the 


order  given  by  (4.10),  which  is  well  within  the  limits  of  a  reasonable  effort. 


The  dressed  coefficients  Oq(4)  and  Y^(4)  are  then  to  be  used  in  the  dressed 
Rayleigh  expansion  (3.5)  to  obtain  the  total  fields. 


5.  Discussion  and  Summary 

The  dressed  Rayleigh  expansion,  i.e.,  (3.1)  with  the  substitution  (3.5), 
takes  the  form  (e.g.,  for  the  electrical  fields  with  the  notation  defined  in 
Fig.  2) 


Eq(x,z)  »  I  { Tq ( i )  pQ  _U)  e 

8,=--“  * 


iCk^x  -  WQ(4)(z  -  g)] 


+  “oU)pO,  +  U)  e 

00 

2.  ( x, z )  -l  Y1  U)  p1  _(4)  e 

oo  * 


i[k^x  +  WQ(  4)  (z  +  g)] 


itk^x  -  W^lMz 


} 

g)3 


(5.1) 


where  YQ(4)  -  6^  e 


iW  (4)g 


E.  ,  and  the  coefficients  a.(4)  and  Y  (4)  satisfy 


Op’  . . * - 0'"'  ~  ’1 

(4.13).  The  convergence  properties  of  (5.1)  can  be  analyzed  for  a  general 
grating  shape  along  the  lines  applied  for  the  SG  in  Section  4.  All  that  is 
necessary  are  the  exact  expression  of  the  corresponding  M  .  and  an  examination 

u>f  X* 

g 

of  its  asymptotic  behavior.  The  explicit  expression  for  M  B  pertaining  to  a 

ul  |  Xi 

16  , 

grating  of  a  general  shape  is  known  (see  also  Appendix  A). 

The  comparison  between  the  convergence  properties  of  the  dressed  and  bare 

Rayleigh  expansion  can  be  explicitly  discussed  in  the  context  of  the  SG.  These 

0  0 

are  determined  by  the  asymptotic  behavior  of  M  0  and  M  .,  respectively.  With  regard 

m,  x  m, x 

to  the  bare  Rayleigh  expansion  (3.1 )»  the  asymptotic  behavior  is  given  by  (4.3) 
and  (4.4).  Since  n(8)  is  monotonlcally  increasing  and  negative  for  8  <<  1  (see 
Fig.  3),  it  follows  that  the  maximum  value  of  8  for  which  M®  .  is  4-convergent 

(Tlf  X» 


(exponentially)  satisfies 


n(B)  -  A  +  B2  +  ln[B/(l  +  A  ♦  B2)]  -  0  , 

or  equivalently,  using  the  notation  of  Petit  and  Cadilhac, 
n  +  1 


(5.2) 


2H 


n  -  1  =  exp[J^-L] 


B 


n2  -  1 


2n 


(5.3) 


The  solution  to  (5.2)  or  (5.3)  is  B  *  0.66  or  n  *  1.86  or  g/d  *  0.1.  This  value 

1 5 

is  in  agreement  with  the  reported  maximum  g/d  ratio  for  which  the  bare 

Rayleigh  expansion  is  found  applicable.  By  contrast,  the  asymptotic  behavior  of 

£,  Eqs.  (it. 7)  ,  yields  convergence  for  all  B!  Note  also  that 

12  7 

alternative  schemes  for  solving  the  light-grating  scattering  problem  ' 

involve  the  bare  coefficients  A q(&)  and  ) .  By  "dressing"  the  coefficients, 

as  articulated  in  (3.6),  the  convergence  properties  of  these  schemes  are 

expected  to  improve.  In  this  regard,  the  dressed  Rayleigh  expansion  provides 

the  simplest  possible  convergent  scheme. 

The  value  n  *  1.86  deduced  from  (5.2)  for  a  SG  is  somewhat  larger  than 
21  2H 

n  -  1 . 5h ,  ’  which  is  the  rigorously  proven  threshold  for  a  perfect  metallic 

SG  (e1  =  -»).  This  difference  highlights  the  distinction  between  singular  and 

non-3ingular  gratings  (i.e.,  smooth  shapes,  finite  c1 ) .  In  the  former  case 

there  are  surface  charges  which  surely  modify  the  boundary  conditions  on  which 

our  discussion  is  based.  We  can  expect  on  physical  grounds  that  when  there  is 

an  infinite  amount  of  charge  at  the  surface  (with  a  {-function  distribution), 

the  point  (x  -  d,  z  -  g)  separating  the  two  adjacent  selvedge  domains  0  S  x  2  d 

and  d  S  x  S  2d  is  singular  (Fig.  2).  Therefore,  two  such  adjacent  selvedge 

domains  become  “disconnected”.  Consequently,  the  Rayleigh  hypothesis  is  no 
,  B 

longer  exact,  the  matrix  M  is  modified,  and  so  is  the  convergence  criterion. 

We  do  not  consider  the  issue  of  singular  fields  any  further. 
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We  comment  now  on  a  related  related  convergence  criterion  for  the  bare 


Rayleigh  expansion,  which  is  presumably  valid  for  any  dielectric  constant  e.) . 


This  criterion  is  based  on  the  application  of  the  steepest  descent  me 


tho<^5 


to 


an  exact  expression  of  AQ ( JL )  or  C^H).  When  applied  to  a  SG,  the  maximum  value 


of  B  for  which  there  is  convergence  satisfies 


24 


3_t.-t  =  en 
n  -  1 


(5.4) 


where  n  is  defined  in  (5.3).  We  understand  the  difference  between  our  result 

(5.3)  and  (5.4)  in  terms  of  the  limited  applicability  of  the  steepest-descent 
25 

method  for  the  SG  case.  Phrased  differently,  the  corrections  to  the  steepest- 

descent  result  for  the  SG  case  are  not  negligible,  and  hence  it  is  not  clear  how 

stringent  condition  (5.4)  really  is.  [The  solution  to  (5.4)  is  n  *  1.54. 

15 

Numerically,  however,  convergent  results  are  obtained  for  higher  values  of  n 
up  to  about  our  result  of  n  «  1.86.]  This  statement  can  be  simply  verified  by 
examining  the  stationary -phase  function  ij)(x)  derivatives  around  the  stationary 
points  of  z  .  We  find  that  although  i|>"(z  )  <  0,  the  fourth  derivative  is 

3  S 

positive  and  has  a  larger  contribution  than  i|»"(z  ). 

s 

In  summary,  we  have  analyzed  the  two  questions  pertaining  to  the  Rayleigh 
expansion:  the  ability  to  continue  it  to  the  grating's  surface  (the  Rayleigh 
hypothesis)  and  its  convergence  properties.  We  have  identified  a  class  of 
gratings,  the  SG  being  one  example,  for  which  the  Rayleigh  hypothesis  is  exact. 

We  have  also  identified  the  cause  for  the  instability  of  the  Rayleigh  expansion 
in  an  actual  calculation  for  gratings  where  the  Rayleigh  hypothesis  is  exact. 

To  remedy  this  instability,  we  propose  a  modified  expansion,  Eqs.  (3.5)-(3.6) 
and  (5.1),  dubbed  as  the  dressed  Rayleigh  expansion.  The  dressed  expansion  has 

presumingly  excellent  convergence  properties  for  deep  gratings.  For  the  case  of  j 

! 

the  SG,  it  is  explicitly  shown  that  the  dressed  expansion  converges  for  an 
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arbitrarily  large  value  of  8  (or  g/d  ratio),  and  the  associated  matrix  to  be 

inverted  is  of  the  order  8irg/d  CEq.  (4.10)].  Our  analysis  is  valid  for  a  frequency- 

dependent  dielectric  function  since  the  whole  discussion  pertains  to  a  fixed 

incoming-wave  frequency  (all  Bragg  scatterings  are  elastic).  The  dressed  Rayleigh 

22 

expansion  is  applied  elsewhere  to  study  the  very  deep  (8  »)  SG  grating  limit. 


Acknowledgments 

We  would  like  to  thank  M.  Hutchinson  and  W.C.  Murphy  for  stimulating 
discussions  and  criticisms.  We  are  grateful  to  A.C.  Beri  for  his  assistance  in 
preparing  the  figures  and  to  R.  Ferron,  J.  Rouse  and  L.W.  Edwards  for  typing  the 
manuscript.  This  research  was  supported  by  the  Office  of  Naval  Research  and  the 
U.S.  Air  Force  Office  of  Scientific  Research  (AFSC),  United  States  Air  Force, 
under  Contract  F49620-86-C-009 .  The  U.  S.  Government  is  authorized  to 
reproduce  and  distribute  reprints  for  governmental  purposes  notwithstanding  any 
copyright  notation  hereon. 


19 


m 


Appendix  A:  Exact  Equations  for  the  (Reflectivity)  kn(l)  and  (Transmitivity) 


C^(&)  Coefficients 


The  exact  infinite  set  of  coupled  \inear  equations  for  the  coefficients  of 


expansion  (3-1)  have  been  derived  by  Toigo,  Marvin,  Celli  and  Hill  for  a  grating 


of  arbitrary  shape.  We  recap  here  their  derivation  for  the  SG  for  the  sake  of 


completeness  and  to  unify  notations. 


We  choose  to  match  the  tangential  components  of  B  and  E.  For  this  purpose 


we  introduce  the  following  triad  of  unit  vectors  which  are  normal  and  tangential 


to  the  SG  grating  profile  given  by  z  =  g  cos  k_x  (the  notation  of  Fig.  2  is  used 

U 


throughout) : 


-»  A 

n  -  -  xB  sin  k  x  -  z 
G 


t  =  -  x  +  zB  sin  k_x 

u 


b  -  s  , 


(A.  1 ) 


where  B  =  k_g,  and  the  hats  denote  the  unit  vectors.  The  normal  vector  n  and 

ll 


the  tangential  vector  t  need  not  be  normalized  for  our  purposes.  By  projecting 


the  field  (3-D  along  t  and  b,  the  following  exact  equations  are  obtained: 


00  ik.x 

l  e  { AQ ( £ )  exp[iW0(H)g  cos  kQx] 


Ci(i)  exp[-iW1  U)g  cos  kQx] }  -  Qt 


(A. 2. a) 


l  e  lA0U)[W0U)  +  Bk^  sin  kQx]  exp[iW0(i)g  cos  kGx] 


( £)  [-W^Jl)  +  Bkj^  sin  kQx]  expOiW^Ug  cos  kQx]}  -  , 


(A.2.b) 


•-  V  v.v.v.w 
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where 

co  ik  X 

Qb  “  -  I  C^iDe  4  exp[-iWQ(i,)g  cos  kQx] 

00  ik  x 

Qt  “  "  I  ^UU-WqU)  +  Bk^  sin  kQx]  e  1  exp[ -iWQ(  4)g  cos  kGx]  . 

(A. 3) 


Equations  (A. 2)  can  be  cast  into  a  set  of  coupled  linear  equations  for  An(X.) 

ik  x  0 

and  C^U)  by  expanding  all  terms  in  the  complete  set  {e  }.  It  is  remarkable 
that  it  is  in  fact  possible  to  further  eliminate  exactly  AQU)  or  C^U),  thus 
breaking  (A. 2)  into  two  (infinite)  sets  of  equations.  This  is  achieved  by  first 
multiplying  (A. 2. a)  by 


-ik  x  +  iWn(m)g  cos  k  x 

Mb  =  t-W0(ra)  +  8km  sin  kQx]  e  (A. 4. a) 

and  multiplying  (A.2.b)  by 

-ik  x  +  iW.(m)g  cos  k.x 

Ht  •  e  "  0  0  ,  (A.H.b, 

adding  the  two  equations  and  integrating  over  the  interval  0£  x  S  d.  This 
manipulation  yields  the  equations  for  C^H).  By  repeating  the  same  procedure 
with 


kfn\  -ikmx  "  iW.(m)g  cos  k„x 

Mb  *  kTTT  C'W1(B1)  "  Bkm  sin  kGx]  6 


(A. 5. a) 


and 


k(1 )  "ikmX  "  lW1(m)g  003  kGX 
k(0)  e 


(A. 5. b) 


an  equation  for  the  A^ ( Jl )  is  obtained.  In  these  manipulations  we  have  used' 
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d  i(ke  -  k  )x  +  iq  cos  k„x 

/  dx[a  +  b  sin  k_x]  e  1 

0  G 

b(k  -  k  )  . 

[a  -  - rr - ]  d  i  V*  (q>  . 


(A. 6) 


where  Jn  denotes  the  Bessel  function  of  order  n. 
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1.  (a)  Schematic  display  of  a  grating  for  which  there  are  four  obvious  domains 


with  constant  dielectric  constant.  The  domains  are  denoted  by  "a",  "b", 

"c"  and  "d".  The  hatched  and  plain  regimes  indicate  the  dielectric  and 
"air",  respectively.  The  maximum  |z |-extention  and 
and  periodicity  of  the  grating  are  denoted  by 

g  and  d,  respectively.  The  Rayleigh  hypothesis  is  exact  for  a  smooth 
grating  of  the  type  depicted  here  (see  Section  2). 

(b)  A  more  complex  grating  for  which  the  Rayleigh  hypothesis  is  exact,  i.e., 
field  expansions  in  domains  "a"  and  "b"  are  identical  (see  Section  2). 

(c)  The  square-well  grating  for  which  the  Rayleigh  expansion  is  not  exact,  namel 
the  field  expansions  in  domains  "a"  and  "b"  are  different. 

2.  The  sinusoidal  grating  (SG),  defining  the  notation  used  in  the  text. 

3.  The  exponential  convergence  function  $(8),  Eq.  (4.8),  pertaining  to  an  SG. 

Note  that  it  is  negative  for  all  6-values. 
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